SARGE: an algorithm for generating QCD-antennas 
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Abstract 

We present an algorithm to generate any number of random massless momenta in phase 



(N 
> 

O 

' space, with a distribution that contains the kinematical pole structure that is typically 

. found in multi-parton QCD-processes. As an application, we calculate the cross-section 

\ of some e+e" partons processes, and compare SARGE's performance with that of the 
uniform-phase space generator RAMBO. 



Considering that many multi-jet processes will occur in future hadron colliders, such as the 
gn! LHC, it is necessary to calculate their cross-sections. A part of the amplitude of these pro- 

r~| , cesses consists of a multi-parton QCD-amplitude, and it is well known ||l| that the leading 

kinematic singularity structure of the squared matrix elements is given by the so-called an- 
\ tenna pole structure (APS). In particular, for n gluons it is given by all permutations in the 

' momenta of 



(1) 



{PlP2){P2P3){P3P4) ■ ■ ■ {Pn-lPn){PnPl) 

where (piPj) denotes the Lorentz invariant scalar product of the gluon momenta pi and pj. 
Actually, it is this kinematical structure that is implemented in algorithms based on the so 
called SPHEL approximation to calculate the amplitudes But it is expected, and observed, 
that the same structure occurs in the exact matrix elements [^, ^. 

For the integration of the differential cross-sections of the processes under consideration, 
the Monte Carlo method is the only option, and a phase space generator is needed. RAMBO 
is a robust and efficient algorithm to generate any number of random massless momenta 
in their center-of-mass frame (CMF) with a given energy. However, RAMBO generates the 
momenta distributed uniformly in phase space, so that a large number of events is needed to 
integrate integrands with the APS to acceptable precision. Especially when the evaluation 
of the integrand is time-consuming, which is the case for the exact matrix elements, this is 
highly inconvenient. 
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In this paper, we introduce SARGE, an algorithm to generate any number of random mass- 
less momenta in their CMF with a given energy, distributed with a density that contains the 
APS. We shall show that it takes account for a substantial reduction in computing time in 
the calculation of cross-sections of multi-parton processes. We briefly sketch the outline of 
the SARGE-algorithm; a fuller discussion, appropriate to hadronic initial states as well, will be 
given elsewhere 0. 

The name SARGE stands for Staggered Antenna Radiation GEnerator, and is inspired by the 
structure of the algorithm. It consists of the repeated use of the basic antenna density for the 
generation of a momentum k, given two momenta pi and p2- 

d^fe.p.;.) = A.XA-'^)«(t»)i Jff^.(i5iiI)Ji^) . (2) 

vr {pik){kp2) \{piP2)J \{P1P2)J 

Here, 5 is a function that serves to regularize the infrared and collinear singularities, as well 
as to ensure normalization over the whole space for k: therefore, g{^) has to vanish sufficiently 
fast for both ^ — > and ^ — > 00. At this point, we take the simplest possible function we 
can think of, that has a sufficiently regularizing behavior. We introduce a positive non-zero 
number and take 

2 log 

which forces the value of ^ to be between (^^^ and and is normalized such that J dA = 1. 
Let us immediately adopt the notation 

^ - ^^^^^ and 6 = M . (4) 



(pm) {P1P2) 

The main motivation to make the regularizing function depend on ^1 and ^2 is that it makes dA 
completely invariant under Lorentz-and scale transformations of the momenta. Consequently, 
the number gives a cut-off for the quotients ^1 and .^2 of the scalar products of the momenta, 
and not for the scalar products themselves. It is, however, possible to relate to the total 
energy ^/s in the CMF and a cut-off sq on the invariant masses, i.e., the requirement that 

{Pi+Pjf>so (5) 

for all pairs of momenta pi 7^ pj. This can be done by choosing 

s _ (n + l)(n-2) 
" ^ 2 ' 

where n is the total number of momenta. With this choice, the invariant masses {pi + k)'^ and 
{k + P2)'^ are regularized, but can still be smaller than sq so that the whole of the demanded 
phase space is covered. The sq can be derived from physical cuts px on the transverse momenta 
and ^0 on the angles between the outgoing momenta: 



So = 2p|. • min ^1 - cos 9q , ^1 + \Jl -p\/s^ ^ 



(7) 
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We now give the algorithm to generate k under the basic antenna density. Let k^, (f) and 6 
denote the absolute value, the polar angel and the azimuthal angle of k in the frame for which 
pi = —p2 with pi along the positive z-axis. To generate k, one should 

Algorithm 1 (BASIC ANTENNA) 

1. determine the direction of pi in the CMF of pi and p2; 

2. generate two numbers ^i, ^2 independently, each from the density g{^)/^; 

3. compute from these the values k^ and cos 9; 

4. generate cf) uniformly in [0, 27r); 

5. construct the momentum k in the CMF of pi and P2; 

6. boost the result to the actual frame in which pi and p2 were given. 

The RAMBD algorithm was developed with the aim to generate the flat phase space dis- 
tribution of n massless momenta as uniformly as possible. The differential density is given 
by 

n 

dVnm) = S{V~s - P')SHP)lld^PiS{p!)9{p^) , (8) 



i=l 



where P = J2i=iPi- denote 



dAl, = dA{qj,q,; q,) , and ^ = " (9) 

To include the APS in the density, one should 

Algorithm 2 (QCD ANTENNA) 

1. generate massless momenta qi and qn in CMF; 

2. generate n — 2 momenta qj by the basic antennas dA\ ,^d^ ,^dA\ ,^ ■ ■ ■ dA, 



n-l . 

n— 2,n' 



3. compute Q = YI!j=i Qj-> the boost and scaling transforms that bring to s/s 
and Q to (0,0,0); 

4. for j = 1, . . . , n, boost and scale the qj accordingly, into the pj. 

This way, the momenta pj are generated with differential density dVn{{p})An^^ {{p}), where 

^QCD.r |^ ^ g(C!,n)g(^l,n)g(^2,n)g('^2,n) " " " g(C-2,'r^)g(C-2,'n) .^qx 

n IWJ 27r"-l ■ {P1P2){P2P3){P3P4) ■ ■ ■ iPn-lPn){PnPl) ' ^ ' 

We point out that, whereas the product dAl ^^ ■ ■ ■ dA'^zl n contains a factor (piPn) in the 
numerator, the scaling transformation carries a Jacobian that is precisely s'^/{piPnf', thus 
leading to a perfectly symmetric APS. 
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Usually, the event generator is used to generate cut phase space. If a generated event does 
not satisfy the physical cuts, it is rejected. In the calculation of the weight coming with an 
event, the only contribution coming from the functions g is, therefore, their normalization. 
In total, this gives a factor 1/(2 log ^m)^""'^ in the density. 

Because we are dealing with gluon momenta, we want to symmetrize the density. This 
can be done by re-labeling the momenta using a random permutation: 

Algorithm 3 (SYMMETRIZATION) 

1. generate a random permutation a £ Sn and put pi <— Po-{i) for all z = 1, . . . ,n. 

An algorithm to generate the random permutations can be found in [^. As a result, the 
differential density becomes 

dK(M) ( ^ E ) ' (11) 

\ perm. / 

where the sum is over all permutations of (1, . . . ,n). An efficient algorithm to calculate a 
sum over permutations can be found in Q. 

When doing calculations with this algorithm on a phase space cut such that {pi +Pj)'^ > sq 
for all i ^ j and some reasonable so > 0, we notice that a very high percentage of the 
generated events does not pass the cuts. An important reason why this happens is that 
the cuts, generated by the choices of g (Eq. dH)) and (Eq- (H)), are implemented only on 
the variables ^^'j that appear explicitly in the generation of the QCD-antenna. Therefore, 
an improvement is obtained as follows. Let Pm denote the subspace of [—1,1]"^ for which 
\xi — Xj\ < 1 for all i, j = 1, . . . ,m, and let us denote the number of .^^'•'^-variables that has to 
be generated = 2n — 4. An improvement is obtained if the generation of these variables is 
replaced by 

Algorithm 4 (IMPROVEMENT) 

1. generate (xi, . . . distributed uniformly in P„^; 

2. define xq = and put, for all i = 2, . . . ,n — 1, 

Because all the variables Xi are distributed uniformly such that \xi — Xj\ < 1, all quotients 
with and {k, I) in {(i — 1, i) , (z, n) | i = 2, . . . , n — 1} are distributed such that they 

satisfy < ^^'■j < ^m- This is an improvement on the previous situation, because then only 
the quotients ^-Zi'^ and ^-i „ with i = 2, . . . ,n — 1 satisfied the relation. In terms of the 
variables Xj, this means that the volume of P^^ is generated, which is + 1, instead of the 
volume of [— Ijlj^^^ which is 2"?. We have to note here that this improvement only makes 
sense because there is a very efficient algorithm to generate the uniform distribution in P^ 



p. The total density changes such that the product of the ^-functions in Eq. (10) has to be 
replaced by 

p jc\^ - 1 ^ Jl if (a^i,... ,x„^) G P„^ , 

- KTlKb^ jo if (x„... ,x.j 0P., , ^''^ 
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where the variables Xi are functions of the variables ^^f^^ as defined by (^2|). Again, only the 
normalization has to be calculated for the weight of an event. 

We compare SARGE with RAMBO in the calculation of the cross-section of the processes 



e ' e 
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qqg, qqqq, qqq'q', qqgg, qqqqg, qqq'q'g, qqggg 



(14) 



The squared matrix element was calculated with the algorithm presented in suitably 
adapted for these processes. We used massless electrons and quarks, and took the sum over 
final-state helicities and the average over initial-state helicities. We also summed over the 
color configurations of the final states. The center-of-mass energy ^/s was fixed to 500 GeV 
for the processes with 5 outgoing momenta, and to 100 GeV for the other processes. The cuts 
on the phase space where fixed with choices of a parameter r, which is related to the cut-off 
So on the squares of the outgoing momenta (Eq. (^)) by 

2ST 



So 



n[n 



I) 



(15) 



where n is the number of outgoing momenta. If r = 1, then sq is larger than the maximal 
value that is kinematically allowed. The couplings and charges in various processes were all 
set to the value 1, since they only contribute a factor to the cross-section, which is irrelevant 
for this analysis. The results of the computer runs are given in the tables below. Presented 
are the final result for the cross-section a in units of GeV~^, the number of generated events 
A^ge, the number of accepted events A^ac; and the cpu-time consumed tcpu in seconds. All 
Monte Carlo runs were performed on a single 440-MHz UltraSPARC-IIi processor, and were 
stopped when an expected error of 3% was reached. 

The final results for the cross-sections are irrelevant in our discussion, and are just printed 
to show that the results with SARGE and RAMBO are compatible within the 3% error estimate. 
The most important conclusion that can be drawn from the results is that SARGE needs less 
accepted events than RAMBO for the given error estimate, especially for small values of r, i.e., 
for phase space that comes close to the singularities of the QCD-amplitudes. (Remember that 
the ratio of the volumes of cut phase space and whole phase space is given by A^ac/Age for 
RAMBO.) As a result, less evaluations of the matrix elements have to be done which accounts 
for a large gain in computer time. It is true that SARGE is "ineffective" in the sense that many 
of the generated events have to be rejected because they do not satisfy the cuts imposed, but 
this is fully compensated by the fact that generating random numbers is much cheaper than 
evaluating matrix elements nowadays. For the last four processes, no results with RAMBO and 
T = 0.01 are presented, but we observe that tcpu > 130, 000 seconds. The fraction of phase 
space covered with five massless momenta and r = 0.01 is 0.893 it 0.001. 



e+e qqg 


T 


0.5 


0.1 


0.05 


0.01 


alg. 


SARGE 


RAMBO 


SARGE 


RAMBO 


SARGE 


RAMBO 


SARGE 


RAMBO 


a 


1.85e-5 


1.85e-5 


1.53e-4 


1.58e-4 


2.61e-4 


2.66e-4 


6.26e-4 


6.41e-4 




7,691 


25, 782 


10, 777 


24,801 


10,806 


37,121 


11,437 


366,614 


A^ac 


5,503 


6,536 


9,436 


20,112 


9,852 


33, 577 


10,860 


359, 447 


^cpu 


251 


293 


429 


899 


451 


1,503 


497 


16,124 
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e'^e — ^ qqqq 


T 


0.5 


0.1 


0.05 


0.01 


alg. 


SARGE 


RAMBO 


SARGE 


RAMBO 


SARGE 


RAMBO 


SARGE 


RAMBO 


a 


9.79e-9 


10.4e-9 


7.72e-7 


7.86e-7 


1.90e-6 


1.83e-6 


7.39e-6 


7.00e-6 




64, 384 


158,678 


32,492 


27, 091 


34, 701 


29, 642 


41,744 


113,368 


iVac 


4,428 


4,551 


9,894 


15,328 


13,081 


22,297 


20, 150 


107,021 




775 


78(3 


L718 


2,(30(3 


2, 25(3 


3, 778 


3, 578 


18, 038 



e+e qqq'q' 


T 


0.5 


0.1 


0.05 


0.01 


alg. 


SARGE 


RAMBO 


SARGE 


RAMBO 


SARGE 


RAMBO 


SARGE 


RAMBO 


a 


5.38e-9 


5.30e-9 


4.07e-7 


4.24e-7 


l.OOe-6 


1.02e-6 


3.95e-6 


3.89e-6 




98, 810 


2 15, 138 


50, 052 


45, 9(33 


(33, 398 


50, 873 


71,254 


36(3, 166 


iVac 


6,696 


7, 022 


15, 392 


25, 883 


23, 989 


38, 145 


34, 584 


345, 323 


^cpu 


1,165 


1,198 


2,664 


4,346 


4,133 


6,434 


5,843 


58, 708 











qqgg 








T 


0.5 


0.1 


0.05 


0.01 


alg. 


SARGE 


RAMBO 


SARGE 


RAMBO 


SARGE 


RAMBO 


SARGE 


a 


1.76e-7 


1.70e-7 


1.86e-5 


1.95e-5 


5.19e-5 


5.27e-5 


5.40e-4 


iVge 


96, 942 


268, 407 


42,321 


86, 608 


50,552 


298, 073 


50,414 


A^ac 


6,579 


7,677 


12,945 


48, 902 


19,091 


223,530 


26,551 


^cpu 


1,363 


1,597 


3,619 


6,398 


3,802 


43,913 


5,287 



c c — qqqqg 


T 


0.5 


0.1 


0.05 


0.01 


alg. 


SARGE 


RAMBO 


SARGE 


RAMBO 


SARGE 


RAMBO 


SARGE 


a 


2.04e-ll 


1.91e-ll 


4.05e-8 


4.08e-8 


1.68e-7 


1.61e-7 


1.48C-6 


iVge 


4,028,648 


4,017, 888 


238, 220 


97, 035 


203,237 


210,325 


176,710 


iVac 


5,616 


5,094 


14,216 


33, 239 


19,522 


121,734 


29, 492 


^cpu 


4, 530 


3,941 


10, 333 


23, 875 


14,159 


87, 756 


21, 407 




e+e qqq'q' g 


T 


0.5 


0.1 


0.05 


0.01 


alg. 


SARGE 


RAMBO 


SARGE 


RAMBO 


SARGE 


RAMBO 


SARGE 


a 


1.05e-ll 


1.05e-ll 


2.19e-8 


2.23e-8 


9.07e-8 


8.86e-8 


7.85e-7 




5, 596, 725 


6,929,475 


436, 225 


188, 693 


377, 384 


522, 602 


305, 426 




7, 730 


8,844 


26, 154 


64, 558 


36, 042 


302, 724 


51,044 


^cpu 


5,882 


6,494 


17, 595 


43, 104 


24, 764 


201,801 


34, 700 
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e+e — 


^ qqggg 








r 


0.5 


0.1 


0.05 


0.01 


alg. 


SARGE 


RAMBO 


SARGE 


RAMBD 


SARGE 


RAMBO 


SARGE 


a 


1.31e-ll 


1.30e-ll 


3.63e-7 


3.54e-7 


1.63e-6 


1.54e-6 


1.85e-5 


iVge 


5,926,016 


6, 650, 538 


366, 538 


131,617 


303, 003 


186,257 


335, 307 




8,194 


8,475 


21,918 


45,157 


29,018 


107,897 


56, 008 


^cpu 


7,407 


7, 398 


18,120 


36,958 


24, 036 


88,318 


46, 673 



As an extra illustration, we also present the convergence to zero of the expected error 
during the Monte Carlo-run for a few cases. In Fig. 0, we plot the relative error as function 
of the number of generated events using a double-log scale. We first of all observe that the 
curves for SARGE are less spiky, which shows that SARGE takes care for a substantial part of 
the singular behavior of the integrand. Every time a RAMBO-event hits a singularity, a term 
much larger than the average so far is added to the Monte Carlo sum, resulting in an increase 
of the expected error. Furthermore, we observe that the SARGE-error converges quicker than 
the RAMBO-error, except in the case of e~^e~ — > qqggg with r = 0.05. However, this is a plot 
of the error as function of the number of generated events, and we know that many SARGE- 
events have to be rejected. A more realistic view is given by a plot of the error as function of 
cpu-time (Fig.^), which clearly shows that SARGE outperforms RAMBO. 
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Figure 1: The expected relative error as function of the number of generated events. 
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Figure 2: The expected relative error as function of cpu-time. 
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